Vortex avalanches and the onset of superfluid turbulence 
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Quantized circulation, absence of Galilean invariance due to a clamped normal component, and 
the vortex mutual friction are the major factors that make superfluid turbulence behave in a way 
different from that in classical fluids. The model is developed for the onset of superfluid turbulence 
that describes the initial avalanche-like multiplication of vortices into a turbulent vortex tangle. 
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Turbulence in classical liquids |l| is governed by the 
interplay of the inertial (the first) and viscous (the sec- 
ond) terms in the r.h.s. of the Navier-Stokes equation: 



dv 

— + V/i = v x w • 

at 



vV 2 v 



(1) 



where u> = V x v, fi = /i + v 2 /2, and /i is the chemical 
potential in the case of an isothermal flow. The transi- 
tion to turbulence is determined by the Reynolds number 
Re = RU I v formed by the characteristic velocity U of the 
flow and the characteristic size R of the system together 
with the kinematic viscosity v. For small Reynolds num- 
bers, the viscous term — uk 2 ^ for a perturbation with a 
wave vector k stabilizes the laminar flow. On the con- 
trary, at large Re ^ 1 the effect of the inertial term 
in Eq. JQ) is dominating, and laminar flow becomes in- 
creasingly unstable towards formation of a chaotic flow of 
eddies. The evolution of turbulence is described by the 
Kolmogorov energy cascade: the kinetic energy of the 
flow is transferred to smaller and smaller length scales, 
via decay into smaller vortex loops along the Richardson 
cascade 2j, until a scale is reached where the energy is 
dissipated by viscosity. 

In superfluids turbulence acquires new features: 
Firstly, because it consists of two inter-penetrating com- 
ponents, i.e., the frictionless superfluid and the viscous 
normal fractions, and secondly, because the vorticity of 
the superfluid component is quantized in terms of the 
elementary circulation quantum k = 2ttH/M where M 
is a mass of a "superfluid particle" (for 3 He, it is the 
mass of a pair of atoms, M — 2m). If both the normal 
and the superfluid components are moving, the turbulent 
state bears more resemblance to the turbulence of clas- 
sical viscous liquids. The superfluid turbulence is now 
a subject of intensive studies, see [3 for a review. A 
new class of turbulent flow becomes possible when the 
normal component is so viscous that it is essentially im- 
mobile while the main flow is supported by the superfluid 
component which contains a large number of quantized 
vortices. We refer to this state as one-component super- 
fluid turbulence. A superfluid motion of such type can 
be realized for flows in narrow channels and in rotating 



superfluids. The velocity of superflow can be character- 
ized by the "superfluid Reynolds number" Re s = U s R/k 
where U s is the mean superfluid velocity with respect to 
the normal component (counterflow velocity). The con- 
dition Re s ~ I corresponds to the Feynman criterion at 
which it becomes energetically favorable to unpin a vor- 
tex from the container wall. If a large nucleation barrier 
for vortex nucleation exists, then vortices are not nec- 
essarily created even at high superfluid Reynolds num- 
bers Re s > I and the superfluid remains in a vortex-free 
counterflow state. For velocities well below the intrinsic 
critical velocity ; 4] v c ~ k/£, superfluid turbulence can 
be initiated if quantized vortices are injected by some 
extrinsic means into the flow. The turbulence develops 
when initial vortices start to multiply and form a vor- 
tex tangle. Recent experiments in rotating 3 He-B 
have revealed qualitatively new conditions for the on- 
set of one-component superfluid turbulence: In contrast 
to the normal fluids, the superfluid turbulence develops 
below certain temperature independent of the rotation 
velocity provided the latter is sufficiently high. Accord- 
ing to the measurements of Ref. |5| , a few seed vortices 
injected into the 3 He B counterflow triggered a transition 
into a state with a turbulent vortex tangle for tempera- 
tures T < 0.6T C while the same injection at T > 0.6T C did 
not create any substantial number of vortices in the final 
state. In these measurements, Re s as high as 200 have 
been reached without any noticeable dependence of the 
transition temperature on the initial counterflow velocity. 
The theoretical arguments supported by numerical simu- 
lations explain Q these conditions in terms of the mutual 
friction between the normal and superfluid components 
that appears in the presence of quantized vortices. 

In this Letter we develop a model that describes the 
onset of one-component superfluid turbulence employing 
the ideas put forward in Ref. H- We consider here the 
initial multiplication of vortex loops immediately after 
injection of seed loops into a superflow with Re s ^ I. 
Since superfluid vorticity is quantized, formation of vor- 
tices during the onset of turbulence is the key issue. We 
consider this process taking into account effects of mu- 
tual friction and derive an equation for evolution of the 
vortex loop density at the initial stage of the transition 
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to turbulence. 

Quantized vortices mediate transfer of momentum be- 
tween the normal and superfluid components which pro- 
duces the mutual friction force on a unit volume of su- 
perfluid 0,0] 

f mf = -a'p s [v s x w 8 ] + ap s [uj s x [u B x v s ]] (2) 

for a normal component stationary in the container 
frame, v n = 0. Here u> s = V x v s and u) s is the unit 
vector. For an array of quantized vortices 



dr n S{r - r„) 



where r„ is the coordinate of the n-th vortex line, and 
5 (r — r„) is a three-dimensional (^-function. The dimen- 
sionless parameters a and a' describe the mutual fric- 
tion between the superfluid and normal components. The 
mutual friction typically arises from interaction between 
quasi-particle excitations and the vortices 0, The 
energy dissipation is determined by the viscous mutual 
friction coefficient a in Eq. (J2J) , and a' is the reactive co- 
efficient. Including the mutual friction force in the Euler 
equation for superfluid velocity we arrive at 

dv 

— - + V/2 = (1 - a')v s x w s + a u> s x (u> s x v s ) (3) 
at 

where a, 1 — a' > 0. 

Regarding the onset of turbulence, the reactive coeffi- 
cient a' in Eq. simply renormalizes the inertial term 
of conventional hydrodynamics [the first term in r.h.s. 
of Eq. |{T{|]. Indeed, performing the stability analysis 
of Eq. we would write v s = v s o + v^., etc., where 
v s o is the solution of a stationary problem while (t) is 
a time-dependent perturbation and obtain a first-order 
in time linear differential equation for v' s (t) with time- 
independent coefficients made out of v s o and its spatial 
derivatives. Making re-scaling of time t = t' /(l — a') 
we arrive at the linearized version of Eq. J2J) where the 
first term has now the prefactor unity like the one in 
the Euler equation [i.e., in Eq. without the viscous 
term]. The coefficient a in front of the dissipative term 
[the second term in the r.h.s. of Eq. J2Jl] is replaced with 
q = a/(l — a') > 0. Therefore, one concludes that the 
first term in the r.h.s. of Eq. drives the flow insta- 
bility towards turbulence in the same way as the inertial 
term in the Euler equation does for a potential flow in 
classical hydrodynamics 0. The fundamental difference 
from classical turbulence is that the dissipative term in 
Eq. (PJ) that stabilizes the flow has now the same scaling 
dependence on velocity and its gradients as the inertial 
term. The relative importance of dissipation in forma- 
tion of turbulence is thus described by the intrinsic di- 
mensionless parameter of the superfluid, q = a/(l — a') 
independent of extrinsic quantities such as U s or R. 



The instability leading to multiplication of seed vor- 
tices can be derived from Eq. ©. Taking curl of the 
both sides of Eq. we obtain the vorticity equation 



du) s 
dt 



(l-a')Vx[v s xu s ]+aVx[ii s x[u s xv s ]] . (4) 



Let I be a characteristic size of vortex loops. In an entan- 
gled vortex-loop state, their three-dimensional density is 
n ~ £~ 3 while the vortex-loop length per unit volume0 



(two-dimensional vortex density) is L = In 



i 2 / 3 . 

Let us express the two terms in the r.h.s. of Eq. (0J 
through the vortex density L. Since u> s ~ kL, the first 
term in the r.h.s. that drives the instability becomes 



(l- Q > s L 3 / 2 ~(l- a ')(Cf s -^ 2 



(5) 



The superfluid velocity in Eq. © is assumed to be 
v s = U s — vq where U s is the counterflow velocity, and 
v ~ k/£ is the self-induced velocity for a vortex loop of 
a length I. The kinetic energy of superfluid grows due 
to the increase in the loop density. The kinetic energy is 
taken from the external source at the length scale R and 
is transferred to smaller scales. The energy flow into the 
unit volume of the fluid is dE/dt ~ ElL + where El is 
the energy of the vortex per unit length. However, the 
energy transfer goes on only during the transient period 
when the vortex-loop density increases. When the loop 
density reaches a value such that U s — vq ~ kL 1 / 2 the 
vortex multiplication saturates and the energy transfer 
stops. If the density happens to become larger, it will 
decrease towards L sat ~ (U s /k) 2 while the kinetic en- 
ergy is returned from smaller to larger scales and, finally, 
back to the external source. In other words, saturation is 
reached when the "turbulent superfluid Reynolds num- 
ber" Re<* urb) = U s £/k is of the order of unity. The condi- 
tion Re s » 1 ensures the separation of scales i R that 
is required for formation of a vortex tangle. The limit 
of large Re s oc is equivalent to vanishing Planck 
constant since n oc h so that the vorticity becomes a 
continuous variable like in classical fluids. However, the 
saturation takes place only in quantum superflviids: L sat 
diverges for classical fluids where k can vanish. This dif- 
ference can be seen in the behavior of the turbulent ve- 
locity in normal fluids 0: the velocity at a length scale 
A decreases as v\ oc A 1 / 3 according to the Kolmogorov 
scaling law. On the contrary, the velocity in a superfluid 
turbulence scales as ve ~ n/£ for smaller I. This shows 
that the "effective circulation" in normal fluid decreases, 
n e ff oc A 4 / 3 . The lower limit on A in classical fluids is 
determined by viscosity at the Kolmogorov dissipation 
scale A ~ i?/R e 1/4 Q. 

The multiplication of vortex loops described by Eq. (JSJ) 
can be understood in terms of vortex collisions and inter- 
connections. Such processes were indeed seen in numeri- 
cal simulations on quantized vortices 0, Hoi HlT | . Recon- 
nection of vortices accompanied by formation of vortex 
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tan gle in normal fluids were considered recently in Refs. 
|l2l 1 1 3| . Each reconnection of quantized vortices takes a 
microscopic time of the order of the quasiparticle collision 
time, much shorter than the times involved in hydrody- 
namic processes. It is accompanied by a small dissipation 
within a volume £ 3 where £ is the core radius of the or- 
der of the supcrfluid coherence length. We consider these 
processes as instantaneous and neglect the corresponding 
dissipation. The rate of increase in the vortex loop den- 
sity should be quadratic in n thus h + — Av r n 2 £ 2 . Here 
v r is the relative velocity of the vortex loops, I 2 is the 
loop cross section, a constant A ~ 1 describes an "effi- 
ciency" of the vortex multiplication due to pair collisions. 
Using definition of L the vortex multiplication rate be- 
comes L + ~ v T l?l 2 . The vortex velocity is determined 
through the mutual friction parameters a and a' such 
that Vi = (1 — a')v s — aCJ s x v s . The relative velocity 
of loops v r is proportional to the longitudinal component 
of Vi,: v r ~ (1 — a')v s which agrees with Eq. 10. 

We assume that vortex multiplication occurs first 
within a certain region in the fluid near the location of 
the seed loops from where the created vortex tangle pen- 
etrates into the rest volume of the fluid. The effect of 
the second (viscous) term in the r.h.s. of Eq. (J1J is to de- 
crease the loop density in the multiplication region by in- 
flating the loops due to counterflow and extracting them 
from this region. The viscous component of the mutual 
friction force leads to variation in the vortex loop length 
I ~ 2ttvl ~ a(U s — vo). The length increases while the 
density decreases if U s — v > 0. When v > U s , the den- 
sity increases since the loops shrink due to the friction. 
Thus the rate of vortex-loop density variation due to the 
viscous component becomes 



-a(U, - v )L 3 / 2 



(6) 



The loops expanding according to Eq. © are extracted 
from the region of their multiplication into the bulk if 
U s > vo or vice versa if U s < vq. This process is accompa- 
nied by dissipation of the kinetic energy. The dissipation 
persists until the equilibrium density is reached. 

As we see, both inertial and viscous mutual friction 
terms, Eqs. (0) and ©, have the same dependence on 
the vortex density, i.e., on the vortex length scale. The 
total variation of the loop density in the multiplication 
region is the sum of the two processes, L = L + + L_. 
Putting vq = k/£ = kL 1 / 2 we obtain 



L = j3 \u s L 3/2 - kL 2 



(7) 



where [3 = A(l - a') - a = (A - q)/(l - a'). 

Equation J7J) looks like the Vinen equation |IJ for su- 
pcrfluid turbulence. However, the difference is that the 
coefficient (5 can now have either positive or negative sign 
depending on the mutual friction parameters. One dis- 
tinguishes two limits. In the viscosity dominated regime 



when < 0, the rate of extraction of vortex loops exceeds 
the rate of multiplication; there is no time for vortices to 
multiply since all the seed and newly created vortices are 
immediately wiped away into the bulk fluid. The num- 
ber of vortices in the final state is essentially equal to the 
number of initial vortices, and the turbulent state is not 
formed. The corresponding stable solution to Eq. Q is 
L — > 0. Equation (JJJ) has another point of attraction, 
L — > oo. This would correspond to a decay of initially 
created vortex tangle in a situation without net counter- 
flow: the vortex loops shrink due to the viscous mutual 
friction force. 

In the low-viscosity regime when f3 > 0, the rate of 
multiplication is faster, and the number of created vortex 
loops is large: Each new vortex loop serves as a source for 
producing more vortices. As a result, an avalanche- like 
multiplication takes place, which leads to formation of a 
turbulent vortex tangle. As the number of vortex loop 
grows, the self-induced velocity increases and finally the 
saturated density L sat is reached. Indeed, Eq. JJJ has 
the stable solution L sat — (U s /k) 2 which corresponds to 
the Vinen equilibrium vortex density for the counterflow 
turbulence 14]. For [3 > 0, the solution to an equation 
of the type of Eq. for a constant U s was found in Ref. 
[lij . At the initial stage, when L <g; L sat , the vortex 
density in the multiplication region varies as 



L/L s 



with the rate 



(L sat /L ) 1/2 - (t/r) 



(8) 



r- 1 = PU 2 /2k . 

In Eq. (JHJ, Lo is the (small) initial vortex density at t = 0. 
The time dependence of L turns to exponential L sat —L oc 
cxp(— t/r) for L — > L sat . The time r decreases with 
lowering the temperature since (3 increases. The energy 
ElL + pumped into the fluid goes first to a reversible 
increase in the kinetic energy of the superflow. A part of 
it is then dissipated due to the mutual friction. 

The overall evolution of the vortex density can be seen 
as an interplay of two processes. The first is the turbu- 
lent instability governed by Eq. 10. The second process 
is expansion of vortex loops from the multiplication re- 
gion into the bulk due to the viscous part of the mutual 
friction force. Equation © taken with the opposite sign 
gives the rate of vortex- loop-density flow into the bulk. In 
this form, it coincides exactly with the Vinen equation 
as derived by Schwarz [Jjj- The solution for a constant 
U s has the form of Eq. © with the characteristic rate 

r- 1 = aU 2 /2 K . 

This time increases with lowering the temperature since 
a decreases. 

The instability develops out of initial seed vortices in- 
troduced into the fluid. In experiments |5| these seed vor- 
tices were produced by the Kelvin-Helmholtz instability 
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at the interface between the A and B phases in rotation 
[l5|| . We note that the inequality q <C 1 is also the con- 
dition for propagation of the Kelvin waves [lil Il7j . Per- 
turbations of the seed vortex lines in the form of Kelvin 
waves can grow and multiply under the action of longi- 
tudinal flows. These processes obviously belong to the 
same class and can also be described within our model. 

Consider a container of a radius R rotating with an 
angular velocity f2. If the counterflow velocity were 
kept constant U s = £IR as in the vortex-free state, the 
density L sat would correspond to the number of vor- 
tices in the container N max ~ TrR 2 L sat ~ N 2 q where 
N eq — 2itR 2 Q/k, is the equilibrium number of vortices. 
Such a big number N max of vortices is unstable, therefore 
it relaxes towards N eq . This decay of turbulence takes 
much longer times than r and is accompanied by polar- 
ization of vortex loops and screening of the counterflow, 
which is beyond the scope of the present paper. In the 
final state, with the vortex cluster formed in the cen- 
ter of the container, the counterflow velocity is decreased 
down to U s ~ k/tq ~ (kQ) 1 / 2 where r$ is the intervortex 
distance for the equilibrium density L eq ~ Q/k. 

In Fermi superfluids and superconductors, the param- 
eter q « (wqt) , where u>q is the spacing between the 
bound states of quasiparticles in the vortex core and t _1 
is their lifetime broadening due to scattering from the 
normal compone nt jal • The ratio a/(l — a') decreases 
exponentially 0, Il8l Il9j at low temperatures. As a re- 
sult, below a certain temperature, an avalanche- like burst 
of vortices in the multiplication region up to the density 
L sat can be achieved, and the well-developed turbulence 
can be realized. However, a large number of vortices 
cannot be created for higher temperatures. In 3 He-B, 
q varies from zero, at T — ► 0, to oo, at T — > Tc. with 
q ~ 1 at T ~ 0.6T C 19]. As pointed out in Ref. @, this 
temperature is in a good agreement with the experimen- 
tal observations of the onset of turbulence. In superfluid 
He-II, the friction coefficient is always small except for 
the immediate vicinity of T\. Thus, a an entangled tur- 
bulent vortex state is usually formed. On the contrary, 
in 3 He-A the friction parameter a is large except for low 
temperatures 0] so that the turbulent state would be 
possible only for very low temperatures. 

To summarize, we develop a model for the onset of su- 
perfluid turbulence taking into account the mutual fric- 
tion between the normal and superfluid components that 
appears in the presence of quantized vortices. The con- 
dition for turbulent instability is velocity independent: 
If the viscous component of the mutual friction force is 
small, which happens at low temperatures, the avalanche- 



like multiplication of seed vortex loops takes place. The 
created vortex loops form a dense vortex tangle which 
further evolves into a developed turbulent flow. On the 
contrary, the number of created vortex loops is low and a 
turbulent state is not formed in the viscosity dominated 
regime realized at higher temperatures. 
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